Grand Unification of the Sterile Neutrino 



Biswajoy Brahmachari and Rabindra N. Mohapatra 
Department of Physics, University of Maryland, College Park, MD 20742, USA 



The simplest way to simultaneously understand all existing indications of neutrino oscillations 
from solar and atmospheric neutrino deficits and the LSND experiment, seems to be to postulate 
a sterile neutrino. We present a realistic grand unified model based on the gauge group 5'O(10) x 
SO(fO)' that leads to the desired masses and mixings for the sterile and the known neutrinos needed 
to understand the above observations while fitting those of the known charged fermions. The model 
is a grand unified realization of the recently proposed idea that the sterile neutrino is the lightest 
neutrino of a mirror sector of the universe which has identical matter and gauge content as the 
"qq ' standard model. The two S'O(f0)'s operate on the two sectors in a mirror symmetric way and are 

, connected by a mixed Higgs representations whose net effect is to connect the superheavy right 

■ handed neutrinos of the two sectors. 
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^ , I. INTRODUCTION 

o\ : 

' A turning point in the search for new physics beyond the standard model may be at hand in the arena of neu- 
, trino physics. The standard model predicts zero neutrino masses but there are now strong indications for neutrino 
^ oscillations (and hence nonzero neutrino masses) from the five solar neutrino experiments (Kamiokande, Homestake, 
Qs^ , Gallex, Sage and Super-Kamiokande 0J^), four atmospheric neutrino observations and the direct laboratory 

CNj ■ observation in the LSND experiment |U Furthermore, to explain all three experiments, three different scales for mass 
" differences (Ato'^ are needed. Since with the three known neutrinos one can get only two independent Atom's, it has 
1^ , been suggested q] that a fourth sterile neutrino be invoked. In the presence of this extra neutrino species ( Vg), one 
can construct several scenarios for solving the three neutrino puzzles ||-||]. In this letter, we will be interested in the 
scheme , where the solar neutrino puzzle is solved via the oscillation of the Ve to Vg using the MSW mechanism 
| pO| and the atmospheric neutrino puzzle is solved via the — Vt oscillation. The solar neutrino puzzle fixes the 
Amg_^ ~ (0.3 — 1.0) X 10~^ eV^, whereas the atmospheric neutrino puzzle implies that Am^_^ ~ 10~^ — 10~^ '^ 
eV^ and the LSND experiment keeps the 0.3 eV"^ < Aml^_^^ < 10 eV^. In this picture the and i^r are nearly 
Oh, degenerate and near maximally mixed. If the existence of the sterile neutrino becomes confirmed say, indirectly by 
^ KARMEN ||ll| observing i'^ — i>e oscillation or directly by SNO neutral current data to come in the early part of the 
next century, a key theoretical challenge will be to construct an underlying theory that embeds the sterile neutrino 
along with the others with appropriate mixing pattern, while naturally explaining its ultralightncss. In this letter, we 
propose a grand unified model which not only explains the smallness of neutrino masses via the seesaw mechanism 
but it also incorporates the sterile neutrino, whose ultralight mass is naturally explained. The model at the same 
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time predicts mass degeneracy between maximally mixed and v^- 

The underlying framework for the present work will be the suggestion that there is a parallel standard model 
^,P|JT^ which is an exact copy of the known standard model (i.e. all matter and all gauge forces identical). The two 
"universes" communicate only via gravity or other forces that are equally weak. At an overall level, such a picture 
emerges quite naturally in superstring theories which lead to Eg, x E'^ gauge theories below the Planck scale with both 
E'g s connected by gravity. Wc hasten to emphasize that despite this apparent promising connection, no vacuum state 
that leads to the details needed for our neutrino model has been discussed to date. In this paper, we will assume the 
sub-Planck GUT group to be a subgroup of E^ x E'^ in the hope of possible future string embedding of our model. 
For alternative theoretical models for the sterile neutrino, see Ref. [ p5[ . 

As suggested in Ref. []l2| , we will assume that the process of spontaneous symmetry breaking introduces asymmetry 
between the two universes e.g. the weak scale v'^f. in the mirror universe is larger than the weak scale v^k = 246 GeV 
in our universe. It was shown in Ref. |l2[ that with this one assumption alone, the gravitationally generated neutrino 
masses [ p^ can provide a resolution of the solar neutrino puzzle (i.e. one parameter generates both the required 
Amg_3 and the mixing angle sin^29e~s — 10^^). 

In this paper, we construct a complete realistic model for known particles and forces and make detailed numerical 
predictions for the neutrino sector within a grand unified scheme that implements the seesaw mechanism. Since the 
simplest GUT model that implements the seesaw mechanism is based on the 5*0(10), we use SUSY S'O(IO) x S'O(IO)' 
as our gauge group with each 5*0(10) operating in one sector. As is well known, the SO(IO) model generally predicts 
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a hierarchical pattern for all fermion masses including neutrinos. However the LSND results in conjunction with the 
atmospheric neutrino data imply that m^^ ~ rrii,^. The challenge is therefore to construct an SO(IO) model which 
will lead to the above non-generic prediction. 

We impose a mirror symmetry between the two SO(IO) sectors so that the field contents as well as the Yukawa 
couplings in the two sectors are identical to each other and all differences between them arise from the process of 
spontaneous symmetry breaking. In order to constrain the model further, we impose the requirement that it conserve 
R-parity automatically without using any extra global symmetries. This ensures that there is a natural cold dark 
matter candidate in the model. We also impose an additional global permutation symmetry S3 which plays a key role 
in ensuring the mass degeneracy between the tau and muon neutrinos. This leads to a three neutrino texture which 
is similar in form to the one discussed in Ref. The connection between the visible and the mirror sector occurs 
via the mixing of the heavy right-handed neutrinos p8[ |. 

The nontrivial nature of the problem arises from the fact that in a GUT framework the neutrino couplings are 
intimately linked to the charged fermion couplings and it is by no means obvious that with a simple set of Higgs 
fields one can make the observed hierarchical pattern of the charged fermion masses and mixings compatible with the 
apparent non-hierarchical mass and mixing pattern for the neutrinos. 



II. THE MODEL: EACH SECTOR 



The fermions of each generation are assigned to the (16,1)0(1,16') representation of the gauge group. We 
denote them by '^e,n,T in the visible sector and by corresponding symbols with a prime in the mirror sector (as we 
do for all fields). The 50(10) symmetry is broken down to the left-right symmetric model by the combination of 
45 © 54 representations in each sector. The SU{2)b, x U{1)b~l gauge symmetry in turn is broken by the 126 126 
representations and we take three such representations (and denote them by Ao,i.2® Ao,i,2)- The role of the these fields 
is two-fold: First, they guarantee automatic R-parity conservation and second, they lead to the see-saw suppression 
for the neutrino masses | p^ . 

The standard model symmetry is then broken by the 10-dim. Higgs fields of which we take three i?o,i,2- As is 
well-known, the 126-dim. representation contains in it left-handed triplets with B — L = 2. Due to the presence of 
the 54-Higgs field S in the model, couplings of type AAS and HHS are allowed and they lead to induced B — L 
breaking vev's vl which give a direct Majorana mass to the neutrinos leading to the so called type H see saw formula 
pO| written symbolically as 
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where vr is the generic vev of the component of 126. As was shown in the detailed minimization of the 
potential in such theories leads to the conclusion that vl is also suppressed by a see saw like formula (i.e. ~ ^tikl 
where A is an unknown parameter in the superpotential). If we choose ~ 10^"* — 10^^ GeV (so that it is not far 
from the GUT scale) and A ~ 0.1 — 0.01, then we get in the eV range. Note that while the second term in Eq. (1) 
arising from the conventional see saw formula leads to a hierarchical mass pattern for the neutrinos, the first term 
has no such obligation. Thus, if we require some neutrinos to be nearly degenerate, the first term has to be given the 
dominant role as we do here. 

A second point is that in the effective MSSM derived from the model, the low energy Higgs doublets will be assumed 
to be linear combinations of the doublets present in all 10 as well as 126 dimensional multiplets. In principle this 
situation can be realized by appropriate arrangement of parameters. 

Next we assume the invariance of the action under a discrete permutation symmetry 6*3 under which the ^, 5't), 
(Ai, A2) and (_ffi, H2) transform as doublets whereas Hq Aq and the rest of the fields transform as singlets. The 
same discrete operates in the mirror sector (i.e. no mirror version of 6*3). This then restricts the form of the Yukawa 
part of the superpotential to the following form: 

Wy = hi^,^,Ho + h2 (*p*^ + ^r^r) Ho + h3 ^,{^^ Hi + H2) + hi [{^^ - *r) i^l 

+ 2 H2] + h *e *e Ao + /2 (^^*^ + ^r^-r) Aq + /3 [(*m*/^ " *r*r) Ai + 2*^^,A2] 

+ /4 *e (*mAi + *r A2) (2) 

Using Eq. 2, we can write down the quark and lepton mass matrices for the visible sector as follows. 
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M^ = Mfo + Mi^26 ; M^'^mS. + M.^ ; M^^ - 3A/i^26 ; = M/j? - SM/'ae (5) 

Even though apriori, it may appear from Eq. 3 and 4 that there are 24 parameters in the visible sector mass matrices, 
the actual number is 16 due to the S3 invariance of the theory which yields eight relations among them. They are 
£3 — £2 ^ 02 ^ aj^ ^ £1 bj^ ^ di _ b2 — dp_ — ^ ^ j^q^ procced to determine the remaining parameters 

cb cs as aa ' 04 C4 ' 63 da 64 d4 be "6 . . . . . 

in such a way they give rise to observed fermion masses and quark mixings and vanishing charged lepton mixings. 
The best fit is obtained for the following values for standard model parameters at the GUT scale (all masses in GeV 
units): nit = 112.00, = 0.370, m^, = 0.0011; nib = 1.115, = 0.0148, to<j = 0.0013; CKM mixing parameters 
si2 = 0.2201; S23 = 0.031; Sis = 0.0039 and the lepton masses nie = 0.00033, = 0.0699, to^ = 1.1817 we find the 
values for the parameters a, 6, c, and d listed in table I. They in turn enable us to determine the Dirac neutrino mass 
matrix for the visible as well as the mirror sector. Although the Dirac mass matrix for the neutrino does not play 
a significant role in the masses and mixings in the individual sector, we will see in the next section that it plays a 
crucial role in the mixing between the two sectors. 



III. CONNECTING THE TWO SECTORS AND PREDICTIONS FOR NEUTRINOS 

The neutrino mass matrix has two contributions as is seen from Eq. 1. For vji near 10^^ GeV, the largest entry 
from the second term in Eq. 1 is of order 10^" eV in the 33 element and much smaller in other places. As far as the 
first term goes, its form is dictated by Eq. 2 and we choose it as follows: 




M,, ^ \ Ai Bi Di ] in eV (6) 

Where ~ This mass matrix looks very similar to the one analyzed in [p^ . In order to obtain the predictions 
for neutrino masses and mixings, we need to know the structure of the mass matrices in the mirror sector and the 
connection between the visible and the mirror sector. 

The exact mirror symmetry between the visible and the mirror sector implies that at the level of the superpotential, 
all couplings in the mirror sector are identical to those in the visible sector. We will assume that the spontaneous 
symmetry breaking breaks the mirror symmetry so that actual mass matrices will exhibit differences. For simplicity, 
we will assume that all doublet vev's in the mirror sector differ by a common ratio from those in the visible sector 
(i.e. v'^f./vwk = C )■ Since this asymmetry will effect the fermion masses in the two MSSM's, we will expect the 
B — L-breaking scales and the GUT scales to be different. This in turn will imply that the induced triplet vev's will 
also be different in the two sectors. Our strategy will therefore be to scale the Dirac mass matrix for the mirror sector 
by a common factor but introduce arbitrary triplet vev's in the mirror sector. 




M,,,, = A„, B„, D„, in eV (7) 



Where, Am — qiAi, A'^^ — qiAi, Dm ~ q2Di and B^. — q2Bi. Let us now try to connect the two sectors which we do 
by postulating the Higgs fields (16, 16') ® (16, 16') (denoted by x®X)- There can now be a connecting term between 
the two sectors given by 

W ^g,^,^',X + 9'A^t.-^'^ + -^rK)x (8) 

We now give a vacuum expectation value to the i/'^v'^' element of x ® X fields. Then only the right handed neutrinos 
of both sectors get connected. This in conjunction with the Dirac masses of both sectors introduces a mixing matrix 
between the two sectors of the following form (assuming for simplicity gc = g'c)' 
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M^^, = g^M^nM-^M^n <x> (9) 

Where, 

/ Ar A', \ 

in GeV (10) 

Where, Ar = hAi, A',. = hA'^., Br = I2B1 and Dr — I2D1. We diagonalize the complete neutrino mass matrix to 
obtain the following absolute values of the mass eigenvalues (in eV's): m^i^ = 90.56, m^'^ = — 90.56, m^/ = 0.0034 

rrii,^ = 1.51, m^^ = -1.509, TOi,_^ = 0.001. The squared mass diffences (in eV^) are Am^.^ = 9.9 x 10^^, Ato^_^ = 
0.003 and Amg_^ = 2.27, where the numbers are given in eV^. The fitted values of the parameters are, a — 0.005, 
Ai ^ 0.0253, A\ = 0.050, Bi = 0.675, A = 1.35 given in eV units, qi = 10, 92 = 60, h = 1.2 10^^, h = 0.5 lO^^ and 
gc < X >— 6.9 X 10^^ given in GeV units. The mixing matrix of the six neutrinos in the basis {ve, v^i, Vt, v'^, v'^, i/'^) 
is approximately given as, 

/ -0.99 0.037 0.039 -0.00025 \ 

-0.031 -0.85 -0.52 -0.00072 
0.019 0.525 -0.85 -0.00043 
-0.042 0.0014 0.00071 -0.999 0.0062 
-0.0053 -0.850 -0.525 

V 0.0032 0.52 -0.85 / 



O" = 



(11) 



Combining this with the mixing angle for the charged leptons, we obtain the final mixing matrix among the four 
neutrinos which looks identical to the corresponding top-left 4x4 submatrix of with only the j/g — i^t entry reduced 
by a factor of 2 because of the presence of a small 13 element in the charged leptonic mass matrix. We have varied the 
parameters of the model to see the allowed range for the rrii,^ relevant for the LSND experiment and find consistent 
solutions for the range 0.5 < nii, jeV < 1.5. Therefore, the j/^ and together could play the role of the hot dark 
matter for the upper allowed range of the masses. Also note that the zeros in the neutrino mixing matrix simply 
means that those entries are less than 10^"*. 

We thus see that in this model not only are all three positive indications of neutrino oscillations are explained but 
the mixing between the heavier sterile neutrinos v'^ and v'^ and the active neotrinos are consistent with all known 



oscillation data such as for example the one from the CHOOZ |21 1. What we find very interesting is that with only six 
parameters describing the entire 6x6 neutrino mass matrix (three active and three sterile) and every other parameter 
fixed by the charged fermion masses, four neutrino masses and 12 mixing parameters that link the active to sterile 
neutrinos which could have observable consequences are all completely consistent with known data. 

Turning now to the consistency of our model with big bang nucleosynthesis (BBN) , we recall that present observa- 
tions of Helium and deuterium abundance can allow for as many as 4.53 neutrino species if the baryon to photon 
ratio is small. There are also a new interesting possibility for generation of Icpton asymmetry in the presence of sterile 
neutrinos which will effect the upper limit on the neutrino number p^ . The relevant parameter that determines if 
extra neutrinos contributes via oscillation of the known ones is Am^sin^26 = Sbbn- It has been argued that in the 
absence of neutrino asymmetries, Sbbn < 10~^ eV^. For large masses of the two mirror neutrinos, this bound is 
not satisfied. We therefore invoke the new mechanism proposed in Ref. to evade these bounds thus restoring 
consistency with the BBN constraints. 

The second feature of mirror universe models is the existence of the mirror photon, which could have experimental 
manifestations. One arena would be the BBN; but as was discussed in Ref. [Q, this problem can be ameliorated 
by the assumption of asymmetric inflation between the two sectors p^ . A second constraint comes from the fact 
that 7 — 7' kinetic mixing p6|| is highly constrained by the BBN considerations. In our model, since it arises from 
the mixed Higgs field x ^^^d thus depends on the splittings among the various sub-multiplets in the mixed field x, at 
the phenomenological level, this will serve to fix the intra-multiplet splitting. It was noted by Carlson and Glashow 
p6| that positronium-mirror-positronium oscillation could also constrain the 7 — 7' mixing. But in our model, we 
assume the mirror weak scale to be about ten times large so that e+'e^' bound state will have mass of about 10 MeV, 
preventing the possibility of this oscillation. 

In conclusion, we have presented a supersymmetric grand unified model for the sterile neutrino which can explain 
the solar, atmospheric and LSND data. The new features of the model are (a) the use of the type II seesaw mechanism 
to explain the smallness of the neutrino masses while at the same time accomodating maximal z/^ — Vr oscillation 
with degenerate neutrinos and (b) the prediction of the neutrino mixings among six light neutrinos (three active and 
three mirror) with very few parameters, which are consistent with all known constraints. 

This work is supported by the National Science Foundation under grant no. PHY-9421385. We wish to thank 
Markus Luty for some discussions. 
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TABLE I. The fitted values of ai, bi, d and di 
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